clear all;
eps = 1.0e-4;
n = 273; %学号32103073
x0=zeros(n,1);
a = unidrnd(2,n,1);
global G b;
G = a * (a') + unidrnd(2) * eye(n);
b = 0.5 * G * ones(n,1);
x_star = FD(x0,eps);     %最速下降
x_star = Newton(x0,eps); %阻尼牛顿法
x_star = BFGS(x0,eps,n); %BFGS方法
x_star = CG_M(x0,eps);   %共轭梯度法
% func   函数子m文件
% gradf1 函数梯度子m文件
% FD     最速下降子m文件
% Newton 阻尼牛顿法子m文件
% BFGS   BFGS子m文件
% CG_M   共轭梯度法子m文件
